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Abstract 

Nontrivial fixed points of the hierarchical renormalization group are computed by nu- 
merically solving a system of quadratic equations for the coupling constants. This ap- 
proach avoids a fine tuning of relevant parameters. We study the eigenvalues of the 
renormalization group transformation, linearized around the non-trivial fixed points. The 
numerical results are compared with e-expansion. 
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1 Introduction 



The renormalization group (RG) is a nonperturbative approach to critical phenomena 
in statistical mechanics and Euclidean field theory. It offers a satisfactory explanation for 
the behavior of statistical models in the critical regime, in particular the appearance of 
universal quanities. It also suggests block spin transformations as an actual scheme for 
computations of critical quantities. Last not least it has proved to be a powerful tool in 
rigorous investigations. Although a lot of efforts have been put into the subject since the 
seminal work of Wilson it is far from being closed. 

The main idea is to think of a critical system in terms of the flow of effective actions 
generated by a block spin transformation. This flow will have fixed points. An instructive 
example is a system with two possible phases separated by a second order transition in the 
absense of an external field. Then we typically find two stable fixed points characterizing 
the different phases and one unstable fixed point signalling a phase transition. Critical 
indices can be calculated as eigenvalues of a linearized transformation at the unstable fixed 
point. Given any critical system the most important task is thus the determination of its 
RG flow fixed points and the flow in their vicinity. 

Let us consider a Euclidean scalar field on a <i-dimensional unit lattice A. The action is 
S(4>) = Sq(4>) + V(4 > ) with V(— <f>) = V{4>). The interaction is called local if the Boltzmann 
factor exp(— V(cfi)) factorizes into a product Y\ xt =\ F((j)(x)). Unfortunately, locality is 
only preserved in an approximate sense by block spin transformations. A main technical 
difficulty therefore is to find an efficient parametrization of the effective actions generated 
by the flow. Also it is not a priori clear what the best choice for the block spin is. A 
guideline is to demand the effective actions to be approximately as local as possible. 

The block spin transformation which we will consider here relies on a decomposition 
of the Gaussian measure d// v (0) determined by the free action So(4>) = \{<p, v -1 </)) into 
a background and a fluctuation part. The free covariance (propagator) is v = (— A) . 
The block spin transformation then consists of an integration step, a rescaling step of 
the background field and the block lattice, and a subsequent subtraction of the field 
independent term. In this setup the non-locality of the effective interaction is traced back 
to the fact that the fluctuation covariance despite its exponential decay couples fields 
further apart than a block distance. 

In the hierarchical model [Q] the free covariance (— A) -1 is replaced by a hierarchical 
counterpart (— A)j^, r . Both covariances share a similar long distance decay, although the 
hierarchical covariance is not translation invariant, and the models therefore behave simi- 
larly in the critical region. As a gain the hierarchical block spin transformation preserves 
locality in a strict sense. A hierarchical RG step turns out to be a nonlinear transformation 
of the local Boltzmann factor F((p), which is a function of a single variable rather than a 
functional of a field configuration. 

One point of view is to think of a hierarchical model as an approximation to its full 
brother. This is obviously only justified if there exists an interpolation between the two 
which admits the computation of systematic corrections around the hierarchical situation. 



Such a scheme is still lacking although a number of first steps into this direction can be 
found in the literature, see, e.g., |J. The main difference between the hierarchical and the 
full model is the absence of wave function renormalization in the former. Wave function 
renormalization comes about in the full model because the effective interaction contains a 
kinetic term as a part of the quadratic piece. This is clearly discarded once the effective 
action is forced to be ultra-local. 

A second point of view is to think of hierarchical models as nontrivial statistical models 
by themselves which are perfectly suited for a RG approach. Due to their nontrivial 
interaction they are far from being exactly solvable. Their phase structure and critical 
behavior is as rich and complicated as in the full case. This point of view seems somewhat 
academic since we do not know of any real ensembles which belong to the universality class 
of a hierarchical model. Real systems do not show this peculiar breaking of translation 
invar iance. 

The third point of view is to consider hierarchical models as a test ground for both con- 
cepts and methods towards a RG treatment of the full model. The problem of non-localities 
in the effective action is believed to be of purely technical nature. This is convincingly 
supported by rigorous studies of many full critical models Q which have been preceeded 
by studies of their hierarchical counterparts. A full RG flow will always contain as a part 
a flow of local quantities which is accurately modelled by a hierarchical flow. 

An exciting possibility within the hierarchical context is to tune the dimension pa- 
rameter d away from integer values. The hierarchical transformation does not involve any 
lattice geometry and contains d and a scale parameter L as variables which can be taken 
real valued. Let us stress that this can be done completely independent of perturbation 
theory. A beautiful picture emerges that shows a sequence of thresholds between two 
and four dimensions where new fixed points appear. For instance, four dimensions is the 
threshold below which an unstable double well fixed point exists which governs three di- 
mensional physics. This picture is believed to be true also in the full model although we 
do not know of a RG formulation which incorporates lattice geometry only in form of a 
dimension parameter which can be tuned at will. 

In this paper we investigate the question of RG flows in the hierarchical model. Usually 
small coupling constants are required to do analytically such an investigation. To study 
infrared fixed points we have no small coupling parameters at hand. Nevertheless, we will 
show that we can do RG calculations for models which are not asymptotically free by 
using some algebraic equations. 

After splitting off the quadratic fixed point, which corresponds to the unordered phase, 
we expand the Boltzmann factor in terms of local operators. We will use two kinds of 
operators, simple powers of fields and normal ordered products. Then the fixed point 
equation becomes a set of algebraic equations for the coupling constants. We truncate 
this system and solve it numerically. This turns out to be much better than searching 
for fixed points by iterating block spin transformations and fine tuning the initial values. 
It allows us to determine even the higher fixed points with multi-dimensional unstable 
manifolds to high accuracy. Then we linearize this set of equations and diagonalize it 
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numerically around the various fixed points. This directly computes the critical indices of 
the hierarchical model without having to rely on any kind of expansions. 

The algebraic equations can also be written down for the full model in terms of a 
polymer expansion. It is less clear then how to truncate this system to a reasonable size. 
But we believe that this is managable and will admit a similar numerical treatment. 

The algebraic equations can be solved for any value of d > 2. Inserting a power 
series ansatz in e for the coupling constants at d = 4 — e, we obtain the e-expansion 
for the double well fixed point and analogously at the other critical dimensions. The 
resulting set of equations can be put into a recursive form which we iterate exactly using 
computer algebra. We compare this e-expansion with the numerical results of the previous 
computations both for the fixed points and the critical indices. 

Finally, we investigate the flow around the fixed points in terms of the eigenoperators 
of the linearized flow. It turns out that the eigenoperators at the double well fixed point 
are approximately normal ordered powers in a small field region but show exponential 
decay at large fields. 

Hierarchical models have been investigated by many authors. Closest to our approach 
is the constructive work of Koch and Wittwer on the nontrivial fixed point ||. The 
bifurcation picture is also described within the framework of a continuous hierarchical RG 
(as opposed to the lattice approach) in the article of Felder ||. Another investigation of 
the e-expansion for the double well fixed point was done by Collet and Eckmann j?]]. 

The present contribution to the subject is an investigation of the flow on the computer 
in terms of an algebraic formulation both for the fixed points, the critical indices, the 
eigenoperators, and the e-expansion. 



2 The Model 

We consider hierarchical RG transformations, defined by 

F^F', F\<P) = |^ 7 (C)F(C + /?# d , 

P = L 1 -^ 2 . (1) 

In Appendix 1, we discuss the hierarchical model as an approximation of a field theoretic 
model and motivate the hierarchical RG transformation. 

d// 7 denotes the 1-dimensional Gaussian measure with covariance 7 > 0, defined 
through 

d/i 7 (x)/(x) = (27r 7 )- 1/2 y_ oo dxexpf- — J f(x). (2) 

L denotes the block size of the corresponding blocking transformation. It is not difficult 
to convince oneself that 

F HT {<t>) = LV(^-i) • exp A (3) 
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is a fixed point of the transformation eq. (|T|). Fut is called high temperature fixed point. 
Let us write the Boltzmannian F product, 



F{4>) = F HT (cP) Z(<f>) . (4) 
The RG transformation of Z is again of the type eq. (||), 

z»z', z' ((/>)= [ dfx Y (c) z(c + /?V) Ld • (5) 



with the changed parameters 7' and (3' given by 

= L - l - d / 2 , i = L" 2 7 . (6) 

In the following, we shall always choose L d = 2. As we shall see below, this choice 
will lead to quadratic fixed point equations. Note, however, that it is not obvious how to 
realize a block spin transformation for such a scale factor on a lattice. For L d = 2, we 
have 

(3 = 2 -( d -2)/( M ) ; = 2 -(d+2)/(2d) _ / 7 ) 

Furthermore, we choose 7 such that j' = |(1 — P' 2 ), i.e., 



7" 1 



:> .2 2 / rf -l). (8) 



This type of conventions agrees with those of ref. J|]. We study in this paper the fixed 
points of the RG transformation eq. (||) that are even in <p, i.e., the even solutions of the 
fixed point equation 

Z{$)= [ d^,(()Z(( + 0cP) Ld . (9) 



Let us remark that if Z is a solution of the fixed point equation (||), then Z$, where 
Zs{4>) := Z(5(f)), 5 > 0, is a solution of eq. (||) if 7' is replaced by j'5~ 2 . The consequence 
of dealing with Z instead of F is that (3 is replaced by (3' < (3. This has important 
consequences in the constructive approach of Koch and Wittwer || . 



3 Computation of Nontrivial RG Fixed Points 

A necessary ingredient for a successful study of the transformation (|l]) is a good choice 
of coordinates. Good here means that already a reasonable number of coupling constants 
provide a good approximation of the full problem. Furthermore, good coordinates should 
be also easy to deal with in practical calculations. 

To study the transformation (|l]), one might consider the expansion of F((p) in powers of 
(j). However, if such an expansion is truncated at a finite order, the resulting approximation 
will probably not define a reasonable (i.e., positive) Boltzmannian. The problem is that 
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the coefficients change their sign from order to order, like in the Taylor expansion of 
exp(— x 2 ). 

A good choice (though perhaps not the best) is to use the split eq. (|j), and consider a 
Taylor expansion of Z((p) in powers of (f>. This choice of coordinates has been successfully 
employed in a rigorous proof that the transformation ([!]) has a nontrivial fixed point in 
d = 3 dimensions || . It has the advantage that the expansion coefficients do not fluctuate 
in sign. Furthermore, the correct large 4> behavior of the Boltzmannian is implemented 
automatically. 

3.1 Expansion of Z in Powers of <p 2 

Let us define a rescaled function P{<p) through 

z(</>) = p (0A/2Y) • (io) 

The RG transformation for P is 

p ^ p', p'(<f>) = J d M i /2 (C) P(C + 04>f • (ii) 



We expand P{4>) in powers of 



™= (12) 



The specific choice of normalization of the pi in eq. ([12]) turned out to be suitable for the 
numerical fixed point solver to be described below. In terms of the expansion coefficients 
pi, the RG reads 

Pl^p'l = J2 S rPmPn, (13) 

m,n 

with 

pmnyl+m+n {2{rn+n))\ ^ Q<l<m + n 

f V8/ (m+n-l) ! sj (21) ! (2m) ! (2n) ! ~ ~ (14) 

else . 

If we look for a RG fixed point, we have to study the infinite set of quadratic equations 

0=pt-J2SrP*mPn- (15) 

m,n 

3.2 Numerical Solution of the Fixed Point Equations 

A straightforward numerical treatment of the problem defined by eq. ( |l5|) becomes possible 
if we truncate the sum over I in eq. (0), introducing a highest index l maX ) 

(max 1 21 

^ = g»?7IBJI- (16) 
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1 


Pi 


h 





0.752806717034 


10+° 


0.423398205312 


10" 


Its 


1 


0.481272697982 


10+° 


-0.113124959171 


10" 


17 


2 


0.313506765870 


10+° 


0.347242613663 


io- 


17 


3 


0.186261032043 


10+° 


-0.179714118083 


io- 


17 


4 


0.100696164171 


10+° 


0.395226286632- 


io- 


17 


5 


0.499270725225 


nr 1 


0.117550777471 


io- 


16 


6 


0.228929876623 


nr 1 


0.570517883570- 


io- 


17 


7 


0.977563729148 


io- 2 


0.622745240266 


io- 


17 


8 


0.390718140134 


io- 2 


0.471372511658- 


io- 


17 


9 


0.146546809430 


10~ 2 


0.307351198867- 


io- 


17 


10 


0.515497714660 


IO" 3 


0.122533903732 


io- 


17 



Table 1: The p\ and ft for I max — 10 in d — 3 dimensions 



Then, of course, also the sums over m and n run only from to Z max , and the fixed point 
problem consists in the study of Z max + 1 quadratic equations, 

fl = 0, 0</<Z max , 

fl = PI- SrPmPn. (17) 

m,n<Z max 

We shall denote solutions of these equations by p\. We used the routine C05NBF from 
the NAGLIB library Q for a numerical solution of the fixed point equations in the range 
2 < d < 4. Notice that the dimension dependence enters through the /^'-dependence of 
the structure coefficients Sf 171 . The program requires an initial guess of the solution. If 
the program is successful, it returns a solution, together with the values of the // for this 
particular solution. In table |I|, we show as an example the pf for Z max = 10 and d = 3, 
together with the /* that can be considered as a measure of the error of the solution. 

We used truncation parameters Z max in the range 10 . . . 50. For the 2-well fixed points, 
I max = 20 was completely sufficient. To get an impression of the finite l max effects, the 
reader is invited to study table 0, where we give the first 11 coefficients of the 2-well fixed 
point in three dimensions for Z max = 10, 20, and 30. 

We define an effective potential V(<f>) through 



v w~-*W=q)- (18) 

In figure 1, we show our results for the 2-well potential for d = 2.1 through d = 3.8 in steps 
of 0.1. In all cases we used l max = 20. The deepest potential corresponds to d = 2.1. With 
increasing dimension, the 2-well gets flatter and flatter until it vanishes in four dimensions. 

From naive power counting and the studies of section ||, one expects that n-well fixed 
points occur when the dimension goes below the threshold d n = 2n/(n — 1). Note that 
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1 


^max — 10 


^max — 20 


^max — 30 





0.752806717034 


10" 


1) 


0.752859732932 


io- 


u 


0.752859732933 


io- 





1 


0.481272697982 


10" 


-0 


0.481191004612 


io- 





0.481191004610 


io- 





2 


0.313506765870 


io- 


-0 


0.313445974082 


io- 





0.313445974081 


io- 





3 


0.186261032043 


io- 


-0 


0.186254920100 


io- 





0.186254920100 


io- 





4 


0.100696164171 


io- 


-0 


0.100729195010 


io- 





0.100729195011 


io- 





5 


0.499270725225 


io- 


-l 


0.499755196476 


io- 


1 


0.499755196485 


io- 


1 


6 


0.228929876623 


io- 


-l 


0.229416956212 


io- 


1 


0.229416956222 


io- 


1 


7 


0.977563729148 


io- 


-2 


0.981866078064 


io- 


2 


0.981866078150 


io- 


2 


8 


0.390718140134 


io- 


-2 


0.394316100256 


io- 


2 


0.394316100333- 


io- 


2 


9 


0.146546809430 


io- 


-2 


0.149410716099 


io- 


2 


0.149410716173 


io- 


2 


10 


0.515497714660 


10" 


-3 


0.536648340900 


10" 


3 


0.536648341727 


10" 


?, 



Table 2: The first 11 coefficients "p\ of the 2- well in d = 3 dimensions for 
three different values of / max 



exactly at these thresholds the operators (j) become relevant with respect to the Gaussian 
fixed point. Thus 2-wells are expected to exist for d < 4, 3-wells for d < 3, 4-wells for 
d < 8/3 and so on. 

By starting our program with a suitable initial guess, we were able to find the 3- well 
potentials. Figure 2 shows the results for d = 2.1 through d = 2.7 in steps of 0.1. The 
deepest potential corresponds to d = 2.1, the flattest one to d = 2.7. In all cases we used 
again Z max = 20. 

The search for the 4-wells was a bit more difficult. It turned out that one needed more 
than 20 couplings for a reasonable parameterization. Figure 3 shows the 4-well potential 
in d = 2.1 and in d = 2.3 dimensions. We used l max = 30 in both cases. 

In all cases, we made a further check on our results as follows: We represented the 
Boltzmannian F by its values on a grid of typically 300 sites on the interval < ^ < A, 
where A was typically 10 or 20. The RG steps were then done by numerically performing 
the integral |], using the NAGLIB routine D01GAF. This routine determines an approxi- 
mation of an integral if the integrand is given on a finite number of points. All fixed points 
determined by the algebraic method were converted to functions on the (/>-grid and then 
checked for stability under iterated application of the integration method. 

4 The Linearized Renormalization Group 

In this section we shall report on a numerical study of the eigenvalues of the hierarchical 
RG, linearized around the nontrivial fixed points. Note that the eigenvalues related to the 
full transformation (|l|) are the same as those related to the transformation (|5|) for Z. The 
reason is that fixed points F* and Z* differ only by the factor Fht that is a fixed point 
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i 


2- well 


3- well 


4- well 





2.1000000 


2.1000000 


2.1000000 


1 


0.4787297 


1.9715290 


1.9851950 


2 


-1.1724335 


0.4611953 


1.0966180 


3 


-3.3445644 


-0.7546772 


0.4105210 


4 


-5.8303542 


-2.2270214 


-5.6565547 


5 


-8.4946711 


-3.9790124 


-1.7706071 



Table 3: The 6 leading eigenvalues a\ at d = 2.1 for the 2- well, the 3- well, 
and the 4- well, respectively 



itself. The transformation for Z around a fixed point parameterized by coordinates pi is 

(Pi + Si) ^ (pf + Si) = J2 ST (P*m + Sm) (Pn + <*n) ■ (19) 



We expand to first order in 5, 



# = (20) 



and identify the linearized RG transformation with the matrix R, 

R ln = 2j2sr n P* m - (21) 

m 

The matrix R is not symmetric. It can, however, be shown that its eigenvalues are real, cf. 



lemma I^J below. We used the NAGLIB procedure F01AKF and F02APF to compute all 
eigenvalues of Aj of R, with i = 0, . . . , Z m ax- It turns out that all eigenvalues are positive 
and can therefore be exponentiated. We define 

Ai = L a * , (22) 

where L denotes the block size, in our case 

Table || shows, as an example, the first five eigenvalues o« of the linearized RG around 
the 2- well, the 3- well, and the 4-well in d = 2.1 dimensions, respectively. Positive eigenval- 
ues a,i are called relevant, negative eigenvalues are called irrelevant, and a zero eigenvalue 
is called marginal. In accordance with the expectation, the n-wells have n — \ (nontrivial) 
relevant eigenvalues. 

One always observes an eigenvalue ao = d that corresponds to the trivial volume 
operator. It corresponds to the fact that the fixed point Z* itself is an eigenvector with 
eigenvalue Ao = 2. The next eigenvalue a\ is related to the critical exponent u, via |j] 

v = — . (23) 
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d 


V 


d 


V 


2.1 


2.08886 


3.0 


0.64957 


2.2 


1.36234 


3.1 


0.62570 


2.3 


1.09916 


3.2 


0.60484 


2.4 


0.95704 


3.3 


0.58640 


2.5 


0.86534 


3.4 


0.56995 


2.6 


0.79985 


3.5 


0.55516 


2.7 


0.74993 


3.6 


0.54182 


2.8 


0.71011 


3.7 


0.52973 


2.9 


0.67729 


3.8 


0.51877 



Table 4: Results for the critical exponent v in the hierarchical model for 
the 2-well 



In table |I| we show our results for the exponent v for d = 2.1 through d = 3.8 in the 
case of a 2-well. It is interesting to compare the d = 3 result with the exponent v for 
the 3-dimensional Ising model. The best known estimates for the latter are in the range 
0.624 . . . 0.630 |i"o|| . For d — > 2, the deviation of the exponent v from the value in the full 
model increases. We know from the exact solution of the 2-dimensional Ising model that 
v = 1 for d = 2. In the hierarchical model we observe v > 2 already at d = 2.1. 

It is interesting to look also at the eigenvectors of the linearized RG transformation. As 
an example we consider the eigenvectors of the transformation eq. (||), linearized around 
the trivial fixed point F*((j)) = 1 and around the nontrivial 2-well fixed point in d = 3 
dimensions. The eigenvectors o[°\(j)) for the linearization around F* = 1 can be given 
exactly: 

OJ O) (0)= const -H 2n (^\ , (24) 
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with 

7 



The Hk denote the Hermite polynomials. We compare these functions for i = 1, 2 with the 
corresponding eigenvectors of the RG linearized around the 2-well fixed point in d = 3, 
see figure 5. The figure shows that the eigenvectors in d = 3 are of similar shape as 
the corresponding functions for the trivial fixed point. Note, however, that the "full" 
eigenvectors contain a factor Fht and thus have a completely different large <f> behavior 
than the Hermite polynomials. 

5 e-Expansion 

In this section, we want to compare the numerical results obtained in the previous section 
with expansions in e = d* — d, where d* is one of the threshold dimensions d n = 2n/(n — 1). 
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This section is organized as follows: In subsection 5.1 a parametrization of the hierarchical 



Boltzmannian is introduced that is suitable for the e-expansion. The recursion relations for 



the expansion coefficients af shall be derived in subsection [5.2|. In subsection p\3| we will 



derive the recursion relations for the computation of the e-expansion for the eigenvalues and 



eigenvectors of the linearized e-expansion. In subsection 5A we shall give some comparisons 
of the e-expansion results with the numerical results. 

5.1 Expansion of F in Wick Monomials 

We consider the expansion of the Boltzmannian F in terms of normal ordered powers 
of cj) 2 . This expansion shall be the basis for the e-expansion to be studied below. One 
expands in terms of normal ordered monomials 



The reason to divide the coefficients o; by j l is that a/ becomes 7'-independent for the 
fixed point. Under a RG step the coordinates o; transform according to 



a, 



i[ = (3 2l Y,Cr n ^ma n . (27) 



The sum in eq. (27) is restricted to |m — n\ < I < m + n, and the 'structure coefficients' 
are given by 

C mn = (2m)! (2n)! ^ 

1 (m + n — l)\ (I + n — m)\ (I + m — n)\ 

Normal ordering with respect to a covariance 7 is defined through 

: exp(a^) :^= exp ^— ^-a 2 + acf^j . (29) 

The normal ordered powers of eft can be expressed in terms of Hermite polynomials, 

7W2 



■~<-\2) MvwJ- (30) 

That we use a direct expansion of F as a basis for the e-expansion seems to be in contra- 
diction to our statement above that it should be much better to use the split F = Fht Z, 
and then expand Z. The e-expansion is, however, an expansion about F = 1 (and not 
about Z = 1). Furthermore, as we shall see from the recursion relations to be derived 
below, the effective Boltzmannians to a given order in e live in a finite dimensional space 
of coupling constants. It is thus irrelevant which coordinates are chosen in this space as 
long as they span this space. 
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5.2 e-expansion for the a/ 

(k) 

In this subsection we derive the recursive relations for the coefficients a ; defined through 



oo 



U = E ^ • ( 31 ) 

k=0 



We determine by e-expansion the coefficients of the infrared fixed point at d = d* — e 
dimensions starting from the trivial fixed point in d* := 21* /(I* — 1) dimensions, where 
I* £ {2, 3, . . .}. The following lemma shows how to compute the af^ recursively. 

Lemma 5.1 Suppose that the coefficients (/? _2 ')( m ) are defined by 

oo 

P~ 21 = 2T* (p- 2l ) {m) e m . (32) 



m=0 



For I* G {2, 3, . . .} and = S^q , we have 



(i) x 2(r 2r ) (1) 

nl*l 



«<5(,(*' a := — 7n*F* ' ( 33 ) 

i* 



and, for all I*, N > 1, 



(TV) = 

2 



X^— E f-2^(/5- 2/ ) (JV - fe) a! fe) +Ecr^- fc) c4 fc )) • (34) 

i* —2 fc = i \ m,n / 



For allN >2 , 

N-l 



(TV) 



2(/3~ 2i *)( 1 ) 



/ _ \ 

- 2[{(3- 2l *) {N) a+ Y,^~ 2l *) {N+1 ~ k)a P) + 

^ k=2 ' 

N-l 

+ ££cre +1 - fe) 4 fe) + 



+ 2 £ cr<e } « 

rrt: m^l* 



(35) 



Proof: The TVth order term of the fixed point equation 

f3- 2l ai = J2cr n a m a n (36) 

is given by 

N N 
2Ti S(/r 2I )( ^- fc ) a (*) = ^ ^ Cra ^) a W . (37) 

fc=0 k=0m,n 



11 



Since = 6i t o and (/3 °)( n ) = for n > 0, we obtain 



2^ />> 



AT-l \ AT-1 

t W + £ (^2,)^-*)^) = 2a iN) + £ £ C J 

k=l I k=l m,n 



r«f fc) «i fc) • (38) 



We have used C™ = 5 m ^. Eq. (|3§|) implies for Z ^ Z* eq. (34). Consider eq. ( p8| ) for the 
case N = 1 : 

2Fa ; (1) =2a[ 1} . (39) 

This implies a£ = a<5^. . To determine the constant a, we have to consider 
the case N = 2,1 = 1* : 



This implies 



2(r 2r ) (1) 4 1) = cf r (a?) s 



eq. (H) for 

(40) 



_ 2(/r 2 '*) (1) 



Let us consider eq. (|38|) for the case N > 3, I = I*: 
N-l N-l 

fc=l fc=l m,n 



■mn (N-k) (k) 



(41) 



(42) 



This implies 



iV-2 

(r 2r ) (JV - 1} « + (r 2r ) (1) a[^ 1} + E (r 2l ') (N - k) aP 

k=2 



N-2 

= 2^craL Ar - 1) a+EE c r« 

m k=2 m,n 



(N-k) (k) 
""n ■ 



(43) 



Thus, 



a 



(N-l) _ 



q-21 

N-2 

E_„ 

mj^l* k=2 m,n 

Using CfSa = 2((3~ 21 *)^ and repl 



i r / ^-2 \ 

N-2 -, 

+ 2 £ C^^a+S^W- 

m: m^/* fc=2 m," J 



+ 



(44) 



acing N by N + 1, we obtain eq. (|35l). □ 
With the following lemma we provide explicit expressions for the e-expansion of f3~ 21 . 
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Lemma 5.2 Suppose that the coefficients ((3 ')( m ) of the e-expansion for the term /3 1 
are defined by 



0-i = 2^ ^{f3~ l )^e m . 



m=0 



Then, we have (J3 ')(°) = 1 and 

m—l 



1 



(g-i)(m) _ 

{P ' ("»-*)! V d* 



d*- m , 



for m > 1 . 



(45) 



(46) 



Proof: We have 



/J-' = 22F2^f( 1 ~-) = 221^ exp I — — I 1 



1 



Expansion of the exponential-function on the rhs of eq. ( |47| ) gives 

1 \ 



0- 



221* 



^ l fl]n2\ n , 



71=1 



n! V 



l ~2, 



Furthermore 



n oo 



A* ) E 



n + — 1 



fc=0 



Insertion of eq. (49) into eq. fl4"S| ) yields 



/r' = 22F 



x 1/ Zln2\ n ^ /n + fc- A ( e ^ n+fc 



n=l 



71! 



d* 



E 

fc=0 



Introducing a new variable m = n + k £ {1, 2, . . .}, we obtain 

oo r m— 1 -i , n „o\ m- k 



i 



m=l L fc=0 

This implies the assertion. □ 



(m — k)\ 



I In 2 

d* 



m—l 
k . 



(47) 



(48) 



(49) 



(50) 



(51) 



5.3 e-Expansion for the Linearized RG 

We shall now consider the eigenvalue problem for the linearized RG equation. The eigen- 
values and eigenvectors are computed by using e-expansion. 

The linearized RG transformation is given by the matrix U(a), with matrix elements 



U nl (a) :=2p 2l J2cr n a r 



(52) 
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Lemma 5.3 Consider the eigenvalue equation for the linearized RG equation 

U{a)b = Xb. (53) 

Suppose that the vector a in d = d* — e dimensions, d* = jiry, I* £ {2, 3, . . .} is given by 
the e-expansion 

a m = £ a^e k , (54) 

k: k>0 

and the e-expansion of U is 

t/ = E u{k) e k • (55) 

k 

Suppose that the e-expansions of b and A are given by 

b= & (m) e m , A= J2 A(n)e "- ( 56 ) 

m: m>0 n: n>0 

Suppose that bm = 5 n: — t and = 2 1 ~~ , for fn G {0, 1, 2, . . .}. Then, we have 

N-l 

A^^EEC"*, (57) 
n=0 m 

and, for I ^ m, 



N-l 



b\ N) = * E E (^- n \m ~ C" n) ) W ■ (58) 

2(2 i* — 2 F J n= g m 

Fwri/iermore, 

^=0. (59) 
The proof of this lemma will be given later. 

Lemma 5.4 f7(a) is symmetric with respect to the canonical scalar product, defined by 

(a,b) := E a " & n • (6°) 

n 

Proof: We have 

(u, U (a)v) = E u n T nl v u (61) 

n.l 

where 

T„, :- DjjH, (62) 

Since T n / = I) n , we have 

(«,Z7(o)u) = (U(a)u,v). □ (63) 

Lemma |5.4| shows that the eigenvalues of U (a) are real. The next lemma shows how to 
compute the e-expansion for the eigenvectors and eigenvalues 



14 



Lemma 5.5 Suppose that the linearized RG group equation is given by the following series 
expansion 

U{a) = U^(a)e n (64) 

n: n>0 

and that the Oth order term is symmetric. Let b be an eigenvector with eigenvalue X 
of U (a), i. e. 

U{a)b = Xb. (65) 
Let be a normalized eigenvector with eigenvalue A*- ** ofU^°'(a), i. e. 

[/(o)( a ) fe (o) = A (o) 6 (o) 5 ( b (o) >6 (o)) = 1; (66) 

where the scalar product ( , ) is defined by eq. (ffify. Suppose that the eigenvalue X^ is 
not degenerate. Then, there exists an e-expansion for the eigenvector b with eigenvalue A 



b = Y,b {m) e m , A = ^A (n V' 



such that 



(6(°),feW) = 0, 

for all N > 0. The coefficients b^ N ' and \( N > are recursively determined by 



N-l 



X (N) = ^( 6 (0) ;[/ (iV-n) b (n), 



n=0 



and 



N-1 



Ln=0 



(67) 
(68) 

(69) 
(70) 



where u 1 - is the component of u perpendicular to M°) . □ 



Proof: The eigenvalue equation implies, for all N, 

(U {m) b {n) - A (m) 6 (n) ) = 0. (71) 

m,n:m+n=N 

Thus, 

N-l 

(C /(0) _ A (o)) 6 (at) = £ (U^-^ - A^-»))6W. (72) 

ra=0 

Since (J7(°) - A(°))6(°) = 0, we may add to U N \ N > on the lhs of eq. © a multiple of 
vector M°) such that eq. ( |6"§| ) holds. Scalar multiplication of eq. (72) with gives 



N-l 

( 6 (o) s (C /(o) _ A (°))bW) = £ (&(°), (tr("-«) - A^-"))6^). (73) 

n=0 
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Since is symmetric, the lhs of eq. ( f73| ) is zero. Therefore eq. ( |69| ) is valid. Since the 
rhs of eq. ( |72] ) is perpendicular to b^°\ eq. (ffO|) can be computed by eq. (72). □ 

The next lemma presents the e-expansion of the linearized RG transformation U (a) . 
Lemma 5.6 Let the e-expansion of (3 21 be 

oo 

[3 2l = 2~wJ2 (P 21 ) (m) e m , (P 21 ) {0) = 1 , (74) 

m=0 

/or Z* G {2, 3, . . .}. Then the e-expansion of U 

U = J2 U (k) e k (75) 

k 

is explicitly given by 

tf=2 1 ^EE C r(^) (M «^ (76) 

r=0 m 

The proof of the foregoing lemma follows immediately from the definitions. 

For the recursive computation of the coefficients for the eigenvalues and eigenvectors 
we need the start values of the recursion relations which are eigenvalues and eigenvectors 
of C/(°). 

Lemma 5.7 Suppose that am = 5 m o- The normalized eigenvectors b^ with eigenvalues 
A (o) fu<f>) 

are 

b^=S m ,rn, A(°) = 2 1 -P, (77) 

for allm£ {0,1,2,...}. 

Proof: We have, using a^} = 5m,o> Cf n = &l,n 

E#> = 2^5^. (78) 

Thus 

(C7(0) 6 (0) )n = 2 i-P 5 2 l-f^ 6 (0). Q (?9) 



We finish this subsection with the proof of |5.3| . 

Proof of Lemma By Lemma |5.5| , eq. (|69|), follows eq. ([[7|) . Since 



U^=2^6 n >, A(°)=2^, (80) 



16 



n 


(n) 

a\ 


n 
a>2 


10 


4.59314 


10 1 


-2.89444 


10 2 


11 


-3.41664 


10 2 


2.26283 


10 3 


12 


2.70284 


10 3 


-1.87364 


10 4 


13 


-2.26138 


10 4 


1.63588 


10 5 


14 


1.99286 


10 5 


-1.50104 


10 6 


15 


-1.84404 


10 6 


1.44371 


10 7 


16 


1.78725 


10 7 


-1.45243 


10 8 



Table 5: Numerical results of a\ and 02 for n from 10 to 16 



we have _ 

(^ {0) -A(°)) ni = 2(2-^-2-^)^. (81) 

Therefore 

(^ (0) -A (0) )^= - A w- Sn*. (82) 
Thus, eq. ( |70| ) of Lemma [5~J implies eq. (p8[). □ 

5.4 Comparison of e-Expansion with Numerical Results 

We evaluated the recursion relations presented in the preceding subsection using the com- 
puter algebra program Maple V Release 2. This allowed us to go to relatively high order. 
We always used programs that computed everything exactly (in the form of analytical 
expressions) and programs that solved the recursion relations numerically. Note however, 
that Maple allows for arbritrary high precision in the numerical computations. It was no 

(n) 

problem to compute the coefficients a) for d* = 4 exactly to sixth order in e. However, 

(n) 

the expressions become quite nasty then. As an example we present the coefficients a, 
up to n = 3 in Appendix 2. The general structure for the expansion around d* = 4 is 
that at a given order n the only nonvanishing coefficients are those with I < 2n. The 
corresponding relation at d* = 3 is I < 3n. The e-expansion is expected to have zero 
convergence radius. However, the series are believed to be Borel-summable. For small e 
even the naively summed low order series can be a quite good approximation. In figure 4 
we show the comparison of the 'true' 2-well potential at d = 3.8 with the 1st and 4th order 
e-expansion. The full line gives the result obtained numerically, and the dashed lines give 
the 1st and 4th order approximations, respectively. 

With the numerical version of the program it was no problem to go to orders like 16. 
In table || we show the intimidating growth of the expansion coefficients when the order 
becomes large. 
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i 


n 


u i,f 





0.5000 


0.5000 


1 


0.0833 


0.0833 


2 


0.0556 


0.0445 


3 


-0.0324 


-0.0190 


4 


0.1468 


0.0888 


5 


-0.5743 


-0.2015 



Table 6: Comparison of the e-expansion coefficients of v in the hierarchical 
model U{ with the ones in the full model v^t 



With the help of the recursion relations of the preceding subsections we also determined 
the e-expansion for the exponent v. We again used an exact version of the program that 
was practicable up to order 6, and a numerical version that could be used to higher order. 

In table we show our results for the expansion coefficients of v, compared with those 



of the full model [11]. The first two orders are exactly the same (the coefficient v\ is 
1/12 in the hierarchical and in the full model). This might be due to the fact that the 
e-expansion of the exponent r\ starts at order e 2 . 

In table ^ we give the results of resummed series for v up to order e k , for k < 5, 
and d = 3.0 . . . 3.8. For the larger values of e, the signal of the divergence of the series 
is obvious. For comparison we also quote our numerical result ('true') and the result of 
a Borel-Pade summation of the 6th order e expansion (BP). The latter was obtained as 
follows: The V{ were divided by i\, and the diagonal Pade approximation of the resulting 
Taylor series was determined. From the resulting rational function Q(e) the estimate for 
v was then obtained by numerically computing 

/•oo 

v(e) = I dt exp(-t) Q(te) . (83) 



There is a quite good agreement with the 'true' results. (For d = 3, the diagonal Pade of 
the Borel transform had a nonintegrable singularity on the positive real axis.) 



6 Conclusions 

In this paper, we have demonstrated that at least for hierarchical models, an algebraic 
computation of fixed points and exponents is feasible. An extension to TV-component 
models and general values of L could be easily done. Of course, many new ideas are 
necessary to do the same thing in full models. The crucial question here is the proper 
choice of parametrization of the Boltzmannian. An interesting question which certainly 
deserves study is whether the e-expansion for the RG flow (in the Wilson sense) in full 
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d 


k = 1 


k = 2 


k = 3 


k = 4 


k = 5 


'true' 


BP 


3.0 


0.6161 


0.6523 


0.6060 


0.8255 


0.3639 


0.64957 




3.1 


0.6007 


0.6293 


0.5961 


0.7228 


0.4235 


0.62570 


0.62599 


3.2 


0.5864 


0.6084 


0.5855 


0.6580 


0.4771 


0.60484 


0.60136 


3.3 


0.5730 


0.5895 


0.5744 


0.6144 


0.5160 


0.58640 


0.58791 


3.4 


0.5605 


0.5723 


0.5630 


0.5837 


0.5369 


0.56995 


0.56973 


3.5 


0.5487 


0.5568 


0.5515 


0.5612 


0.5423 


0.55516 


0.55525 


3.6 


0.5377 


0.5428 


0.5402 


0.5440 


0.5379 


0.54181 


0.54183 


3.7 


0.5274 


0.5302 


0.5291 


0.5303 


0.5289 


0.52973 


0.52973 


3.8 


0.5177 


0.5190 


0.5186 


0.5189 


0.5187 


0.51877 


0.51877 



Table 7: Results of resummed series for v up to order e fc , for k < 5. 'true' is 
the numerical result, BP is the result of a Borel-Pade summation explained 
in the text 



models could be performed to an order that is competitive with what has been done in 
the conventional framework. 
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Appendix 1: Field Theory and Hierarchical Models 

This appendix discusses the hierarchical model as an approximation of a field theoretic 
model and motivates the hierarchical RG equation. 

The generating functional for Greens functions of a scalar field theory on a ti-dimensi- 
onal continuum R rf is given by the following (formal) infinite-dimensional integral 

Z[J\ := N f J] <W(*) exp -A0)] exp [-V(<f>) + (J, <f>)} . (84) 

is a scalar field, and J is an external source on H d . V{4>) = J ze ^d V(4>(z)) is the 
local interaction term, and \{4>, — A</>) is the free part of the action. A is the Laplacian, 
and N is a normalization factor chosen such that 

dnytf) := AT I] exp [-±(0, -A0)] (85) 

zeR d 
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is a normalized Gaussian measure, v is called propagator. The canonical scalar product 
is defined by 



<f>(z)i/>(z), 



(86) 



for fields (f>, ip on R d . 

For RG calculations it is more convenient to use generating functionals with external 
fields ip, 

Z{^) := [ d[x v (<t>) exp [-V{4> + V)] • (87) 



The two generating functionals are related by 

Zty) = Z[J] exp [-|(J,vJ)] 



For the definition of the hierarchical model, we introduce the notion of a hierarchical 
lattice (or multigrid). For L G {2,3,...} and j G Z divide R rf into hypercubes of side 
length dj := L~ J a, where a is a unit length. Denote the set of all these hypercubes by 
Aj. Aj can be considered as a lattice with lattice spacing a,- by identifying the centers of 
the hypercubes with the lattice sites. The location of the hypercubes can be chosen in 
the following way. For a hypercube y E Aj let y be the open hypercube of y. For y E Aj 
and x E Afc we suppose that y n x = or y C x. In the latter case we write y£x. The 
hierarchical approximation is given by the following replacement 

v — ► J2 v ^ 



where, for z, z' G R d , 

yj(z,z') :=a?- d 7 ^- (90) 

Here, x and x' are the uniquely determined hypercubes of Aj such that z G ar, z' G x'. v J is 
called fluctuation propagator. The above replacement represents the fact that in general v 
can be decomposed into a sum of fluctuation propagators which are exponentially decaying 
with decay length aj. This exponential decay is simulated by the Kronecker delta on the 
rhs of eq. (p0|). There are other ways to define the hierarchical approximation. But all 
hierarchical approximations share the property that the fluctuation propagators v 3 have 
compact support. Insertion of replacement (^) into eq. (87) yields, using the convolution 
formula of Gaussian measures, 



n i 



iP) exp 



(91) 



We define an ultraviolet cutoff by setting the propagator v J = if j > n and an infrared 
cutoff by setting \ k = if k < j : 



,. n i- n 



i=j+i 



i=j+i 



(92) 



20 



The effective generating functions and zj™ \ are related by 



An), 



(93) 



Since we started with a local interaction V, the effective generating functions obey the 
following factorization property 



zpw= n 4 n w). 



(94) 



Since the kernel of v 3 is constant on hypercubes of Aj, we can assume that the fields tp 

are constant arguments of Zj , i.e., do not depend on z G R d . Therefore, eq. (|93| ) is 
equivalent to 



n 

y.ye_x 



dii vj {4>)z\ n) {y\4> + ^) 



(95) 



for x G Aj_x- For translation invariant models Zj n \y\(f)) does not depend on y G Aj. Let 
us define 



Zf\4>) :=Z 



(96) 



From eq. (|95|) and definition eq. (96) follows, using ip — > a,_iip, 



d H {^)zf\^ + L l -H) 



(97) 



For F' := [zf\) L ~ d and F := {zf ] ) L ~' , we obtain the hierarchical RG transformation 
eq- (|). 

Appendix 2: Some Results of e-Expansion for the af 1 



R 
T 

l o 

l 

l 
,( 2 ) 



In 2 
1 

864 

12 R 2 (3 - 2) - R 2 (19 - 18) T 



-31104T + 41472 



i? 2 



(98) 



216T-432 
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-2R 2 (7R -8) + R 2 (7R- 12) T 
17280T- 24192 

R 

~Tu 

12 R (3 R - 2) - R (19 R - 18) T 

-10368 T + 13824 
R (154 R 2 + 1161 g - 459) - 12 g (8 R 2 + 69 g - 27) T 
- 746496 T + 1057536 

R 2 

2592 T — 2592 

-8i? 2 (i?- 1) + i? 2 (5i?-6)T 
72576 T - 103680 

_R2_ 
41472 

-12 R 2 (3 - 2) + i? 2 (35 i? - 18) T 
-1492992 T + 1990656 
R 3 

373248 T - 373248 

R* 
17915904 

(99) 
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Figure Captions 



Fig. 1: Our results for the 2-well potential for d = 2.1 through d = 3.8 in steps of 0.1. 
The deepest potential corresponds to d = 2.1. With increasing dimension, the 2-well gets 
natter and flatter until it vanishes in 4 dimension. In all cases we used Z max = 20. 

Fig. 2: The fixed point 3-wells for d = 2.1 through d = 2.7 in steps of 0.1. The deepest 
potential corresponds to d = 2.1, the flattest one to d = 2.7. In all cases we used l max = 20. 

Fig. 3: The fixed point 4- well potentials for d = 2.1 and d = 2.3. The potential with the 
deeper wells corresponds to d = 2.1. In both cases we used Z max = 30. 

Fig. 4: Comparison of the 'true' 2-well potential at d = 3.8 with the 1st and 4th order 
e-expansion. The full line gives the result obtained numerically, and the dashed lines give 
the 1st and 4th order approximations, respectively. 

Fig. 5: Comparison of the two leading eigenvectors of the RG transformation linearized 
around the trivial fixed point and the 2-well fixed point in d = 3. The single well shaped 
functions correspond to the relevant eigenvalue. The double well shaped funtions cor- 
respond to the first irrelevant eigenvalue. The full lines belong to the non-trivial fixed 
point. 
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